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Abstract 

In selective cooperative relaying only a single relay out of the set of available relays is activated, hence 
the available power and bandwidth resources are efficiently utilized. However, implementing selective 
cooperative relaying in time- varying channels may cause frequent relay switchings that deteriorate the 
overall performance. In this paper, we study the rate at which a relay switching occurs in selective 
cooperative relaying applications in time-varying fading channels. In particular, we derive closed-form 
expressions for the relay switching rate (measured in Hz) for opportunistic relaying (OR) and distributed 
switch and stay combining (DSSC). Additionally, expressions for the average relay activation time for 
both of the considered schemes are also provided, reflecting the average time that a selected relay remains 
active until a switching occurs. Numerical results manifest that DSSC yields considerably lower relay 
switching rates than OR, along with larger average relay activation times, rendering it a better candidate 
for implementation of relay selection in fast fading environments. 
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I. Introduction 

Cooperative relaying has been recentiy proposed as a means of achieving the beneficial effects of 
diversity in wireless communications systems, without employing multiple antennas at neither the receiver 
nor the transmitter. Its operation is based upon the concept of employing wireless relaying terminals that 
assist the communication between a source and a destination terminal by receiving the message sent from 
the source, and then processing it appropriately and forwarding it to the destination. The relaying terminals 
may be either fixed, infrastructure-based terminals placed at selected spots in urban environments aiming 
at extending the coverage while avoiding the infrastructure cost that the deployment of base stations 
entails, or mobile, hand-held devices. In any of the above cases, apart from the apparent robustness 
against small-scale fading, cooperative relaying offers resihence against large attenuations due to path- 
loss, as well as shadowing. This, together with the advantages that cooperative relaying offers in the 
various levels of the open system interconnection (OSI) protocol stack, renders the cooperative concept 
a strong candidate for utilization in future wireless networks [1]. 

The most common cooperative relaying protocols were introduced in [2], where the term "cooperative 
diversity" was used so as to emphasize the diversity advantages of relay employment. In the same 
work, an outage analysis of the cooperative diversity concept was also conducted, showing remarkable 
performance benefits as compared to the case without relaying. Nonetheless, the analysis in [2] concerns 
the scenario where a single relay is available for cooperation. In cases where multiple relaying terminals 
are utilized, the error performance can be dramatically improved if the multiple relay transmissions occur 
in orthogonal channels and are combined by a maximal ratio combiner (MRC) at the destination. However, 
orthogonal channel utilization results in a reduced overall spectral efficiency. To this end, activating only 
a single relay out of the set of available relays has been shown to be an effective means of achieving 
cooperative diversity while hmiting the negative effects of orthogonal relay transmissions, offering thus 
a good tradeoff between error performance and spectral efficiency. 

Previous works on single relay selection in cooperative relaying scenarios include [3]-[6], where the 
relay selection was based on a maximum signal-to-noise-ratio (SNR) policy, thus attaining a diversity 
order equal to the number of available relays, i.e., the same diversity order as for the case where all the 
relays are activated, yet with considerably higher spectral efficiency. Due to the somewhat opportunistic 
usage of the available resources, the relay selection protocol based upon the maximum SNR rule is termed 
opportunistic relaying (OR). A simpler alternative to OR, for the case of two available relaying terminals, 
is the so-called distributed switch and stay combining (DSSC) protocol proposed in [7]. According to 



May 12, 2009 



DRAFT 



2 



DSSC, a single relay remains active for as long as the corresponding SNR is greater than a predetermined 
threshold value; should this condition be violated, a relay switching occurs. The DSSC protocol hence 
requires only a single end-to-end channel estimation for each transmission, reducing thus the overall 
complexity while achieving the same outage performance as OR, albeit inferior error performance [7]. 
In the sequel, we use the term "selective cooperative relaying" to refer to the OR and DSSC protocols, 
i.e., to protocols where single relay selection takes place. 

A. Motivation 

Despite the above benefits of selective cooperative relaying, however, a major issue that needs to be 
addressed is the rate at which a switching of the active relaying terminal occurs in practical scenarios 
where the fading in each of the links involved is time-varying. In fact, this rate reflects the number of 
times per second the system has to switch from one relaying terminal to another, and corresponds to 
a complexity measure regarding the implementation of selective cooperative relaying in practice. More 
specifically, frequent relay switchings may cause synchronization problems due to the fact that the system 
needs to repeat the iiutiaUzation process each time the active relay changes, in order to re-adapt to the 
chaimel conditions of the new branch. Apparently, such synchronization readjustment leads to increased 
implementation complexity, as well as potential delays and outages which may cause severe information 
loss with ultimate deteriorating effects on performance. 

Particularly for the case of DSSC, where only one end-to-end branch is estimated in each transmission 
period, relay switchings have a negative impact on channel estimation along with synchronization. This 
is because each time a relay switching occurs, a previously idle relay needs to be "awakened". Hence, a 
new training sequence needs to be initiated, which may not be long enough to provide accurate chaimel 
estimation, resulting in detection errors in addition to those owing to weak chaimel conditions. To the 
best of the authors' knowledge, the concept of relay switchings in time- varying fading environments has 
not been addressed in the Uterature. 

B. Contribution 

In this paper, we study the effect of the time-varying nature of fading channels in selective cooperative 
relaying applications. In particular, we provide closed-form expressions for the relay switching rate 
(measured in Hz) of OR and DSSC, as a function of the average channel gains and the maximum 
Doppler frequency of each of the source-relay and relay-destination links involved. These expressions 
consider the case of independent but not necessarily identically distributed (i.n.i.d.) Rayleigh fading 
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channels, and account for both AF and decode-and-forward (DF) relaying. Particularly for OR, we derive 
the relay switching rate for arbitrary numbers of participating relaying terminals when operating over 
independent and identically distributed (i.i.d.) Rayleigh fading channels, showing that this rate is an 
increasing function in the number of relays. In addition to the switching rate, we obtain closed-form 
expressions for the average relay activation time, which is defined as the average time interval that the 
selected relay remains activated until the system switches to another relay. A set of numerical examples 
is provided showing that DSSC results in a considerably lower relay switching rate, as well as a larger 
average relay activation time than OR. These results indicate that DSSC may be preferable in fast fading 
scenarios where the channel gains change rapidly making OR difficult to implement due to the frequent 
relay switchings, despite DSSC's inferiority in terms of error performance [7]. 



The rest of this paper is organized as follows. The mode of operation of OR and DSSC, together with 
some basic definitions are given in Section JI] Section |lll] provides expressions for the relay switching 
rate and the average relay activation time of OR, when operating over Rayleigh fading channels, while 
the corresponding expressions for DSSC are presented in Section|IVl Some numerical examples are given 
in Section |Vl Finally, Section |Vl] concludes the paper. 



We consider the cooperative relaying setup where a source terminal, S, communicates with a destination 
terminal, D, with the aid of L relaying terminals which are denoted here by Ri, i G {!,..., L}. The 
relaying terminals may operate in either the DF or the AF mode. In the former case, the relays fully 
decode the received signal and forward a noise-free symbol to the destination, while in the latter case the 
relaying terminals are used as simple analog repeaters which amplify the received signal and forward it 
to the destination without demodulating it. Additionally, the relays are assumed to be half-duplex, in the 
sense that they cannot receive and transmit simultaneously; instead, the source-relay and relay-destination 
transmissions are assumed to occur in time-orthogonal channels. 

In this paper, we adopt the general notation, aAB, to denote the channel gain of the link between 
terminals A and B, so that the channel gain of, e.g., the S-Ri link, is represented by asR^- The fading in 
each of the links involved is assumed to be Rayleigh distributed, with probabiUty density function (PDF) 
given by 



C. Outline 



II. Background 




(1) 
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where Qab represents the average squared channel gain of the A-B Unk, e.g., ^SR, = E [ogj^] with 
E [•] denoting expectation. We denote the maximum Doppler frequency of the A-B Unk by J^ab, e.g., the 
maximum Doppler frequency of the S-Ri Unk is denoted by J-'sRi ■ Moreover, all terminals are assumed 
to transmit with identical power, denoted by Pt, while the noise power in all of the Unks involved is 
identical and denoted by Nq. 

Throughout this work, two selective cooperative relaying protocols are considered: The opportuiustic 
relaying (OR) protocol presented in [4], and a variant of the two-relay distributed switch and stay 
combining (DSSC) protocol proposed in [7], which is referred to as DSSC-B here; this notation is 
adopted since DSSC-B varies from DSSC in exactly the same way as SSC-B varies from SSC-A in [8]. 
The reasoning behind studying DSSC-B instead of the original DSSC protocol presented in [7] lies in 
tUe fact that DSSC-B yields less frequent relay switchings than DSSC, similarly as SSC-B yields less 
frequent switchings than SSC-A [8]. The modes of operation of OR and DSSC-B are given in detail in 
the ensuing subsection. It is worth mentioning that in each case, only a single relay out of the set of 
available relays is activated and denoted by Rf,, resulting in a somewhat distributed version of selection 
combiiung for OR; and a distributed version of SSC-B for DSSC-B. The selection is performed at the 
destination terminal, which collects the chaimel state information (CSI) of all the Unks involved. Then, 
after determiiung the "best" relay, the destination sends a feedback message to the relays indicating the 
activation of the selected relay and the deactivation of the previously selected relay. All other terminals 
remain inactive until they receive a proper activation message from the destination. 

A. Mode of Operation of Schemes Under Consideration 

1) Opportunistic Relaying (OR): The OR protocol consists of selecting a single relay out of the set 
of L available relays, particularly the relay with the highest of some appropriately defined metric, which 
accounts for both the S-Ri and Ri-D Unks and corresponds to performance measures of the ith end- 
to-end path. Under the assumption of equal power transmitted by the potential relaying terminals, such 
metrics may be the min equivalent defined as 

ai = min {asR^ ,aR^D) (2) 
or the "half harmonic mean" equivalent defined as 

Qi = ; . (3) 

asRi + aRiD 

The min equivalent metric accounts for determiiung the end-to-end path based on the weakest intermediate 
link. It is appropriate for DF relaying because it corresponds to an outage-equivalent of the end-to-end DF 
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channel, since tlie outage probability of DF relaying equals the cumulative distribution function (CDF) 
of ai evaluated at the outage threshold SNR. The half harmonic mean equivalent is appropriate for AF 
relaying since it corresponds to a tight approximation of the end-to-end SNR of the S-Ri-D link [9]. 
It should be noted, however, that both criteria lead to approximately the same results, since the min 
equivalent represents a tight upper upper bound of the half harmonic mean equivalent, particularly when 
the SNRs of the source-relay and relay-destination links are very different, e.g., asR^ 3> qr^d [10], [11]. 
For this reason, in the sequel we adopt the max-min criterion for determining the selected relay, Rf,, so 
that the relay associated with the maximum "bottleneck" of the source-relay and relay-destination Unks 
is selected, i.e.. 



2) DSSC-B: The DSSC protocol [7] appUes to the case where there are two relays available for 
cooperation. Its simpUcity over OR lies in the fact that in each training period, only a single end-to-end 
channel has to be estimated. This estimation is used so as to check whether the active branch is of 
sufficient quality. The system thus switches from one relay to another only if the equivalent SNR of the 

active branch lies below a predefined switching threshold, T. Using the min equivalent metric defined 
above, a relay switching occurs when Ofe < -^/T/T, where F = Pt/Nq denotes the ratio of the power 
transmitted by the source and each relay divided by the noise power, i.e., the common SNR without 
fading. 

Nevertheless, there may exist transmission periods where both the available end-to-end channels are 
not strong enough; in such cases the system switches continuously from one relay to another, though 
without reaching the desired SNR level. Therefore, in order to avoid these excessive, as well as unavailing 
switchings, throughout this paper we study a variant of the DSSC protocol denoted here by DSSC-B. 
According to DSSC-B, a switching occurs whenever the SNR of the active branch down-crosses the 
switching threshold, T; the system then stays connected with the new branch, regardless of whether the 
SNR of the new branch is greater or lower than T, until this SNR down-crosses T. Mathematically 
speaking, denoting by a] the min equivalent of the ith branch for a transmission period j, the active 
relay, RI^^, for the ensuing transmission period is determined by 



b = arg max min {asR^ , aR-o) ■ 
ie{i,...,L} 



(4) 




(5) 



< y^T/T or 4 > Vr/r 



If Ri = R2, Ri 



R2 if 



= < 



Ri if 



4~^ > VtJ^ and 4 < Vrjf 
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It is interesting to note that, under tlie min equivalent criterion, the end-to-end path between source and 
destination in both OR and DSSC-B is treated as a virtual channel with channel gain = min {asR, , clr^d)' 
The OR scheme can be thus regarded as a virtual selection diversity scheme, where the instantaneous 
chaimel gain of the ith input branch is a^; the DSSC-B scheme can be interpreted as a virtual SSC-B 
scheme with chaimel gains ai and 02. 

B. Basic Definitions 

The rest of the paper focuses on deriving expressions, as well as providing numerical examples, for 
the following performance measures for both OR and DSSC-B: 

• The relay switching rate, defined as the number of times per second that a switching of the active 
relay takes place; that is, the destination stops receiving from a certain relaying terminal and coimects 
with another one. 

• The average relaying activation time, defined as the average time duration that the selected relay 
remains activated, starting from the time it receives an activation message until the system switches 
to another relaying terminal. 

III. Relay Switching Rates and Average Relay Activation Time of OR 

A. Two Relays, i.n.i.d. Fading 

Let us first consider the two-relay OR scenario, where the fading in all the intermediate links involved 
is i.n.i.d. 

Theorem 1: The relay switching rate of OR with two available relays and i.n.i.d. fading is given by 



SRoR 



(Ol + ^2?'^ {^SR, + ^R,d) {^SR, + ^R,d) 

(6) 



+ 



'K^/W^Q^ ^R2D^SRi^^SR2^SR2 + ^RiD^R^D + ^RiD^R^D^ ^SRi^sr^ + ^SR2^l 



•2 

SR2 



{Ql + n2f^ i^SRr + ^RrO) {^SR2 + ^R2d) 

where = E [a|] = ^SRi^RiO/ {^SRi + ^Rio) represents the average squared value of a^, i G {1, 2}. 
Proof: Let us consider the random processes 

ai (t) = min {asR, (t) , gr^d (t)) (7) 

02 {t) = min {asR2 (i) , gr^d (*)) (8) 
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which correspond to the fading processes of the virtual end-to-end channels S-Ri-D and S-R2-D, 
respectively. Additionally, let us define the random process Z (t) as 

Z (i) = ai (t) - a2 (t) (9) 

so that the active relay in each time instance t is determined by the signum of Z {t) in this time instance, 
i.e., Ri is active at time t if Z (t) > 0; R2 is active at time t if Z (t) < 0. Consequently, in order to derive 
the average relay switching rate it suffices to evaluate the average number of times the process Z (t) 
crosses zero. This is equivalent to obtaining the level crossing rate (LCR) of Z {t), evaluated at zero. 
Then, the positive-going LCR of Z {t) corresponds to the average number of times the system switches 
from i?2 to Ri , while the negative-going LCR of Z (t) accounts for the average number of times the 
system switches from Ri to i?2- 

The relay switching rate equals the sum of positive-going and negative-going zero-crossing rates of 
Z (t), which can be expressed as [12] 



SRoR 







f (0,z]dz+ zf(0,z]dz (10) 



where / (^^^^^ denotes the joint PDF of Z (t) and the time-derivative of Z (t), Z [t). Because of the 
independence of the fading process and its time derivative, in each of the intermediate links involved, 
owing to the Rayleigh fading assumption, the processes Z (t) and Z (t) are independent [13]. Therefore, 



/ \^z, zj can be expressed as / yz, zj = fz (z) f ^ \^)' '^^^^ fz (") and f ^\^] denoting the PDFs of 
Z (t) and Z (t), respectively. Equation (fTOl ) thus yields 

/■O 

SRoR = fz (0) ' 



—00 
00 



\^zj dz+ zf^(z]dz 



(11a) 

= 2/z(0) zf^[z)dz (lib) 

where we used the fact that the two integrals in (II lab are equal to each other since, apparently, the 
number of times the system switches from Ri to R2 equals that of switching from R2 to i2i, in the long 
run. The PDF of Z [t) evaluated at the origin is derived as (see Appendix A) 



(1^1 + ^2?/^ 
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Moreover, the second term in ([TT]) is derived as (see Appendix B) 







"4 [') '^^ = 



i^SRi + ^RiD) {^SR2 + ^R2d) 



^R2D^SRi\J^SR2^SR2 + ^RiD^R^D + ^RiD^R2D\/^SRi^SRi + ^SR2^SR2 



(13) 

The relay switching rate of OR for the case of two relays with i.n.i.d. fading channels is obtained by 
substituting (fT2l) and (fT3l) into (II Ibb . completing the proof. ■ 
Corollary 1: The average relay activation time for OR with two available relays and i.n.i.d. fading is 
given by 

where SRqr is given in 

Proof: Let us denote by pf^ the steady-state probability of selecting relay Ri in the OR setup, i.e., 
pOR _ Pj.|q^ > Q,2}; = Pr{ai < 02}. Then, given that the average switching rate from Ri to R2 
equals the average switching rate from R2 to Ri, the average relay activation time is derived as 

Using the fact that pf^ = 0,i/ {Qi + 0,2) (see eqs. (gOll, dS)), the proof is complete. ■ 

B. Two Relays, i.i.d. fading 

Corollary 2: The relay switching rate of OR with two available relays and i.i.d. fading is given by 

TT 



'^^OR 



•^SRi + •^SR2 + y-^SRi + -^la-D + y-^RiD + ^SR2 + V -^RiD + ^R2D 



(16) 



4^/2 

Particularly for the case where the maximum Doppler frequencies are also identical (i.e., J^sR, = ^R,D = 
T, i ^ {1, 2}), the relay switching rate is given by 

SR^i)\ = ttT. (17) 
Proof: The proof follows directly from Q after simple algebraic manipulations. ■ 
One may note that the relay switching rate in the i.i.d. case is independent of the channel amplitude 
and is determined only by the maximum Doppler frequency in each of the links involved. Moreover, it 
is interesting to note that ([TT] ) yields a relay switching rate which is identical to that of conventional 
selection diversity (i.e., where no relaying takes place) with identical Rayleigh fading, as given in [13], 
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except for a factor of \/2. This implies that the OR setup can be considered as a distributed selection 
diversity scheme with two virtual Rayleigh channels, where the maximum Doppler frequency of each 
virtual channel equals \/2-tiines the maximum Doppler frequency of the intermediate links. This is due to 
the fact that in our case we deal with four time-varying links, whereas in conventional selection diversity 
there are only two time-varying links involved. 

Corollary 3: The average relay activation time of two-relay OR for the i.i.d. scenario with identical 
maximum Doppler frequencies is given by 

ATSii = ^- (18) 
Proof: The proof follows from (fT4l) . in conjunction with ([TT] ). ■ 

C. L Relays, i.i.d. Fading 

Let us now consider the versatile case of L available relaying terminals. For simplicity of the analysis, 
it is assumed that all the S-Ri and Ri-D links experience i.i.d. fading, as well as identical maximum 
Doppler frequency E . 

Theorem 2: The relay switching rate of OR with L available relays, i.i.d. Rayleigh fading and identical 



maximum Doppler frequency in each of the intermediate channels involved, is given by 

SR^n = V2L{L- 1) 7 ^) {j^) ' ^9) 

Proof: Similar to the case of two available end-to-end branches, the switching rate for the L-branch 
case is evaluated through the use of the process 

Z (t) = ai (t) - ak (t) (20) 

where (t) is the fading process of any of the virtual end-to-end channels as defined in ^ and 

Ofc (t) = max a, (t) . (21) 
je{i,...,L} 

It is important to note that, due to symmetry, the rate at which the system switches from Ri to any 
other relay is not affected by the index i and equals half of the overall relay switching rate, since all 
relays are selected with identical probability. Consequently, the statistics of Z (t) are the same for each 
i G {1, L}, hence the relay switching rate for this case is obtained as 

SR^i = Lfz (0) r zf^ (z) dz. (22) 



'We note that the results can be extended to the case of i.n.i.d. fading and identical maximum Doppler frequency. However, 
the resulting expressions would be too complicated and are out of the scope of this paper. 
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Using trivial integrations and the expressions for fz (0) and f ^ yzj given in Appendix C, in ( |49l ) and 



(l52l ). respectively, (l22l ) yields ( fT9l ); the proof is thus complete. As a cross-check, one may notice that for 
L = 2, (O reduces to ([T7]). ■ 

It is interesting to note from ( [T9l ) that the relay switching rate is an increasing function of L, implying 
that, as expected, the larger the number of relays the more frequent relay switchings occurs. 

Corollary 4: The average relay activation time for the i.i.d. scenario with arbitrary number of available 
relays is derived as 

^TSi = (23) 
Proof: Using the same approach as in Corollary [B it follows that 



OR Q jyiid Q rjiid 



where we used the fact that, due to symmetry, pf^ = 1/L. 



IV. Relay Switching Rates and Average Relay Activation Time of DSSC-B 
Theorem 3: The relay switching rate for the two-relay DSSC-B setup is derived as 



T(ni+n2) / _2T_ 



J^DSSC [-1 ) = — 



T 2T T 2T T(ni + fi2) T(fii+2n2) T(2rii+n2) 

2grni _ grni _j_2er'"2 — e'""^ — 2e rriiS"!2 -|- g rniii2 -|- g rniQ2 



T 2T T 2T T(ni+Q2) T(fii+2n2) T(2Qi+n2) (^4) 

2grni _ grni _j_2grn2 _ grna — 2e rriin2 +g rC!iQ2 -)_ g rniQ2 _2 



T 2T T 2T T(ni+ri2) T(fii+2n2) T(2fii+n2) 

2grni _ groi _j_ 2g rn2 — grn2 — 2e rr!ir!2 +g rsiiQa -|- g rnjna —2 



Proof: In the two-relay DSSC-B scheme, a relay switching occurs whenever the negative-going slope 
of the output SNR crosses the switching threshold, T; equivalently, a relay switching occurs whenever 
the process Oj (t), i e {1,2}, crosses \/T/T in a negative-going direction; recall that T denotes the 
common SNR without fading. The switching rate for relay Ri, i G {1,2}, is given by 

(T) = r \x\ 4^ (Vt/T, x) dx (25) 

where ^ (•, •) denotes the joint PDF of a {t) and ai (t). Using (l34l) . (|42)) and the fact that a (t) and 
aj (t) are independent processes, (1251 ) yields 
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The steady-state probability of activating Ri in the DSSC-B scheme is derived in Appendix D as 



[T/T) [F„. (T/r) - 1] 


(r/r) + F2 (T/r)" 




F,. (T/r) [F,. {T/T) - 1] 

Fal (T/r) [ 


1 - (T/r) 

F,. {T/T) - 1] 


+ 2i^i^ (T/r) 
"i - F,. (T/r) 


+ i^ai {T/T) 

+ F2 (T/r)^ 


[F,. (T/r) - 1] 

(27) 


(T/r) [F,. (T/r) - 1] 


1 - 2F,. (T/r) + 2F2 (T/r)" 


+ F,.(r/r) [F,. (T/n-i] 



(28) 



where 

{x) = l- exp (-x/fli) (29) 
denotes the CDF of of {t), i G {1,2}. Therefore, the relay switching rate can be expressed as 

SRdssc {T) = pf (T) Ni (T) + (T) iVs (F) . (30) 



Substituting (|26l)-([29l) into (|30l ) and after some algebraic manipulations, we arrive at (|24] ). ■ 
Corollary 5: The average relay activation time for DSSC-B is given by 

Proof: The proof is identical to that of Corollary [T] We note that the factor 2 in the numera- 
tor is present because the system switches to Ri from the other available relay at a rate that equals 
SRdssc (F) /2 since, apparently, the number of times the system switches from R2 to Ri equals the 
number of switches from Ri to R2, in the long run. ■ 

V. Numerical Examples and Discussion 

A. Implementation Issues 

In classical diversity communication systems (e.g. conventional selection combining or switch and stay 
combining, where no relaying takes place), the switching between the diversity branches causes several 
problems, such as "an internal outage" due to the corruption of the receiver filters and data signal chains, 
as well as phase estimation failures [13]. In cooperative relaying systems, however, where signals have 
to be exchanged between spatially distributed nodes, an additional important issue has to be addressed; 
such issue is time synchronization, which requires the local clocks of the relay nodes to be synchronized, 
requiring various degrees of precision [14]. 

It should be pointed out that in classical switched diversity systems, keeping time synchronization 
between the transmitter and the receiver after switching to a different diversity branch is not a serious 
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problem, since the relative time delays of the different branches are practically identical, because the 
link distances are practically identical. On the contrary, in cooperative relaying systems, especially in 
those employing mobile relays, the time delay between the transmitter and the destination may change 
dramatically as soon as a relay switching occurs. Time synchronization becomes particularly challenging 
in cooperative networks since the network dynamics such as propagation time or physical channel access 
time are in general non-deterministic, because of the relative distances between source, relays, and 
destination. In other words, switching to another relay would provoke a time synchronization readjustment 
between the transmitting relay and the destination, as well as the source. Therefore, it is evident that 
frequent relay switchings are not desirable as they significantly increase the implementation complexity 
of relay selection protocols, rendering them difficult to implement in fast fading scenarios. 

B. Numerical Results 

In this subsection, we illustrate some numerical results, regarding the relaying switching rate and 
average relaying activation time of both OR and DSSC-B. All the results were also confirmed by 
simulations. As already mentioned, the relay switching rate and the average relaying activation time 
constitute a vital part in the design and implementation of selective cooperative relaying, due to the 
complexity issues that each relay switching entails. Fig. [T] depicts the normalized switching rates of both 
OR and DSSC-B versus the ratio TsrJ J^r^d for L = 2 relays, when the fading powers in the S-R 
and R-D links are unbalanced. Specifically, the average channel gains of the S-R\ and S-R2 links are 
assumed equal to each other, and the same is true for the R\-D and R2-D links; the difference in the 
corresponding average SNRs is lOdB, so that V^sRi = ^^RiD i lOdB, i G {1,2}. As observed from 
Fig. [U the switching rate of DSSC-B is always smaller than that of OR, with the latter to be under 
specific conditions even ten times greater than the former. This verifies the intuition that DSSC-B leads 
to less frequent relay switchings, accounting thus for simpler practical implementations. It should be 
noted that in Fig. [T] the switching threshold, T, used is that threshold which leads to the maximum 
possible switching rate of DSSC-B, which is determined through numerical optimization methods. In 
other words, the DSSC-B curves depict the worst case in terms of switching rate, demonstrating that 
even in this scenario the DSSC-B switching rate is much less than that of OR. 

The effect of the number of available relays on the switching rate of OR for the i.i.d. case with equal 
maximum Doppler frequencies in all of the links involved, is plotted in Fig. |2] As expected, the switching 
rate continuously increases as the number of the relays increases, yet in a non-linear manner. The main 
result extracted from Fig. |2] is that in cases where the number of available relays is large, leaving some 
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of them out of the selection set may be preferable in practical applications with high Doppler spread, 
despite the apparent performance cost. 

The normalized average relay activation time of the relays for the same conditions as in Fig. [1] is plotted 
in Fig. |3] Again the average activation time in DSSC-B is the minimum possible, i.e., the parameter T 
is properly selected so as to account for the worst case. Fig. [3] manifests that the average relay activation 
time of DSSC-B is considerably larger than that of OR, regardless of the assumptions on the relative 
S-R and R-D channel strengths. The effect of the number of relays on the activation time for OR is 
shown in Fig. |4l One may observe that the activation time reduces in a non-hnear fashion as the number 
of relays reduces. 

Finally, in Fig. |5l the switching rates of OR and DSSC-B are plotted for L = 2 relays and several 
values of the switching threshold, T, assuming that the channel gains in the S-R and R-D links are 
unbalanced, as well as that TsR^ = '^^RiD- For each of the scenarios shown here, the switching rate 
of DSSC-B is considerably lower than that of OR. The ratio of switching rates is, however, crucially 
dependent on the switching threshold of DSSC-B, implying that the difference of the switching rates of 
OR and DSSC-B shown in Fig. [T] is significantly expanded for T values different from the worst (in 
terms of relay switching rate) case for DSSC-B. 

VI. Conclusion 

We conducted a study of selective cooperative relaying in time-varying fading channels. In particular, 
we derived the relay switching rate for selective cooperative relaying, which reflects the concept of how 
frequently a cooperative scheme switches from one relay to another. Together with the relay switching 
rate, a closed-form expression for the average relay activation time was also derived, which corresponds 
to the average time which a relay remains activated until a switching occurs. Numerical results indicated 
that for the case where there are two relays available, selecting the active relay in a max SNR-based 
fashion (a.k.a. opportunistic relaying - OR) results in a considerably higher switching rate than selecting 
the relay in a DSSC-B fashion. Therefore, considering the complexity issues associated with frequent 
relay switching, DSSC-B leads to a simpler implementation of selective cooperative relaying than OR, 
and may be preferable in fast fading scenarios, despite being inferior in terms of error performance. 
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Appendix A 
Evaluation of fz (0) 

The PDF of Z (t) evaluated at zero reflects the probability that the absolute difference of the virtual 
fading gains of the S-Ri-D and S-R2-D channels lies in the infinitesimal interval [0,(iz]. Therefore, 
fz (0) can be expressed as 

POO 

fz (0) = / fa, (x) fa, (x) dx (32) 

Jo 

where fa, (•) and fa, (•) denote the PDFs of the processes oi (t) and a2 (t), respectively. Using ([T]) and 
(|2]), the CDF of Oj (t), i e {1,2}, is expressed as 

pOO /"OO 

Fa^ (x) = 1 - / fsR, (w) duj 1 fR^D {i^) du 

J X J X 

, / ^SR, + ^R,D 2\ ,,,, 

= 1 - exp X (33) 

V ^''SrMr.d J 

while the PDF of ai{t) is derived by differentiating (|33] ). yielding 

^^SR, + ^RiD ( ^SR, + ^RiD 2 



fa, (x) = 2x — ■ ---- exp . (34) 

Note that the process (t) is also Rayleigh distributed; its average squared value is denoted by Qi and 
is given by 

^SR, + ^R,D 

Then, fz (0) is derived as the integral of the product of two Rayleigh distributions, yielding 



Using [15, eq. (3.321.5)], dMl) yields ([Ell. 



Appendix B 



Derivation of / (•) 

z 

The time derivative of Z {t), Z (t), equals the difference of the time derivatives oi (t) and 02 (t). 



Before proceeding in deriving \^zj, we first obtain the PDF of Oj (t), i G {1, 2}, as 

A (x) = Fr {asR, < ai?,D} (a;) + Pr {a^i?, > clr^d} f„ {x) (37) 

where (•) and (•) denote the PDF of the time derivatives of asR, and qr^d, respectively. 

Given that qsr^ and ur-d are Rayleigh distributed, (•) and (•) can be expressed as [12] 

fasn = ~ "^""P f ^ I ^^^^ 
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/■ i^) = exp ( %- 



(39) 



i.e., asR^ and aji^n are zero-mean Gaussian random variables (RVs) with standard deviations ctsr^ 



T^^SRi yJ^SR, and ctj?^d = T^^R.Dy/^hijD^ respectively. The steady-state probabilities Pr {asR^ < aR^o} 
and FiiasR, > or^d} are given by 



poo 

Plc{asR, < Ur^d} = / fasR- (x) [1 - i^an o (x)] dx 

Jo 

Fr:{asR, > aR^o} 



^SRi + ^RiD 



SR, 

^SR^ + ^R,D 



(40) 
(41) 



The PDF of Oj (t) is therefore derived as 



\/^SR,^SRi 



+ 



\/^RiD^RiD 



(42) 



where i G {li2}. Consequently, it follows from Q that the PDF of the process Z {t) is expressed as 
[16, ch. 6] 

ix) = / r {y + x) r (y) dy. (43) 



Substituting (l42l ) in ( |43l ) and making use of the integral [15, eq. (3.323.2)] 

(~) 



CO / 

exp exp 



ay = . TT exp 



Ol + 02 



Ol + 02 



(44) 



we derive the PDF of Z (t), / (•), as 

z 



Z ^7r3/2 (ytsR, + ^R,d) i^SR, + ^r,d) 

^RiD^SR2 exp 



f^5fli^^Si?2 exp 



^R.D^R^D + ^R^D^R^D 



+ 



^R2D^SRi exp 



+ 



'2n^{QsR2^ln2+<^'^ion.,n) 



+ 



'^SR^^SRr + ^R2D^R,D 

0,R^D^R^Dex.p -irrm r^^To — T 

^SR^J'sR, + ^SR2^SR2 



^SR2^SR2 + ^RiD^RiD 



(45) 



Having an expression for / (•), the last term in ([TT]) is derived using [15, eq. (3.321.4)] as shown in 

z 



(SI. 
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Appendix C 

Relay Switching Rate of OR with L available relays 

Evaluation of fz (O); Let us first derive the PDF of ak {t) as the PDF of the maximum of L — 1 i.i.d. 
Rayleigh RVs, yielding 

L-l L-1 

fa. (^) = E (^) n (^) = - 1) (^) i^t~^ (46) 

;=i j=i 

with fa (x) = (2x/fi) exp [-x^/^) and (x) = 1 - exp [-x^/^l), where J] = Q-skJ"^ = ^r^d/"^ = 
... = = Q.ji^i:,/2 denotes the average squared channel gain in each of the virtual end-to-end 

channels involved. Then, fz (0) is derived as 

poo poo 
fz (0) = / fa (x) fa, (X) dx = {L-l) fl (x) [Fa (x)]^"' dx. (47) 

Jo Jo 
Using the product expansion 



1 — exp 



x' 



( l?7b yields 



Derivation of f {■): Because of the i.i.d. assumption, the PDF of the time-derivative ai (t) is derived 

z 

from (|42]l as 

that is, ttj is a zero-mean Gaussian RV with standard deviation = 'KT\f^. The PDF of the time- 
derivative of a/c (t), Cfc (t), is derived as 

which implies that, due to symmetry, the PDF of the time-derivative of the maximum of L — 1 i.i.d. RVs, 
equals the PDF of each of the L - 1 RVs. Using dH, the PDF of Z {t) = Ui {t) - (t) is derived as 

1 f \ 
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Appendix D 

Steady-State Relay Activation Probabilities of DSSC-B 

Considering that relay switchings in DSSC-B are determined in exactly the same way as branch 
switchings in SSC-B, the steady-state relay activation probabilites for DSSC-B are derived through the 
Markov states of SSC-B given in [8]. Specifically, the Markov chain of DSSC-B yields six states, as 
follows. State 1 corresponds to the case where is active and the overal SNR down-crosses T"; state 
2 corresponds to is active and the overall SNR is below T"; state 3 con^esponds to is active and 
the overall SNR is greater than T"; states 4, 5, 6 refer to the case where R2 is active, and are defined 
analogous to 1, 2, 3, respectively. The stationary probabilities, Tj-j^-S'^^C-B^ j ^ {l^...j6}, of the above 
Markov states are taken from [8, eq. (21)], yielding 

^DSSC-B ^ ^ ^^.^^ '^^'^2 .,9,9 (53a) 

vr?^"^-"^ = :--^r . . . (53b) 

vr.f ^"^-^ = T ... (53c) 

(53d) 
(53e) 

^uaau-B ^ -ill-ill^ ■izj'riz .^2f) 

' [qi + 92] [1 + 2qiq2] - [qi + q2? - 2qlql 

where qi = (T/F); q2 = Fg^ (T/F). Considering that Ri is active in the states 1, 2, 3, while R2 
in states 4, 5, 6, the steady-state relay activation probabilities of DSSC-B are derived as p^^^^ = 
E|=ivrf p^''^ = E?=4vrf yielding m and m. 



B _ 


[1 - qi] qi [1 


- 92] q2 






[qi + 52] [1 + 2qiq2] - 


[qi + q2f 


- 2qlql 


B _ 


Qiq2 [1 - 


-q2] 






[qi + q2] [1 + 2qiq2] - 


[qi + q2? 


- 2qlql 


B _ 


[1 - qif [1 


- 92] 92 






[qi + ^2] [1 + 291^2] - 


[qi + q2? 


- 2qlql 


B _ 


[1 - qi] qi [1 


- 92] qi 






[qi + q2] [1 + 2qiq2] - 


[qi + q2? 


- 2qlql 


B _ 


qhi [1 - 


-qi] 






[qi + q2] [1 + 2qiq2] - 


[qi + q2f 


- 2qlql 


B _ 


[1 - qi] qi [1 


- 92] q2 






[qi + 52] [1 + 2qiq2] - 


[qi + q2? 


- 2qlql 
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Fig. 1. Relay switching rates of OR and DSSC 
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Fig. 2. Relay switching rate of OR versus the number of available relays 
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Fig. 3. Average activation time of OR and DSSC 
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Fig. 4. Average activation time of OR versus the number of available relays 
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Fig. 5. Relay switching rates of OR and DSSC versus the normalized switching threshold 
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